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Statistics of pressure and pressure gradient 
in homogeneous isotropic turbulence 

By T. Gotoh 1 AND R. S. Rogallo 2 


The statistics of pressure and pressure gradient in stationary isotropic turbulence 
are measured within direct numerical simulations at low to moderate Reynolds 
numbers. It is found that the one-point pdf of the pressure is highly skewed and 
that the pdf of the pressure gradient is of stretched exponential form. The power 
spectrum of the pressure P(k) is found to be larger than the corresponding spectrum 
Pc(k) computed from a Gaussian velocity field having the same energy spectrum 
as that of the DNS field. The ratio P(fc)/Pa(fc), a measure of the pressure-field 
intermit tence, grows with wavenumber and Reynolds number as —R^ 2 log(k/kd) 
for k < kd/2 where kd is the Kolmogorov wavenumber. The Lagrangian correlations 
of pressure gradient and velocity are compared and the Lagrangian time scale of 
the pressure gradient is observed to be much shorter than that of the velocity. 


1. Introduction 

The pressure field plays an important role in the turbulent motion of an incom- 
pressible fluid. The pressure, or more precisely its gradient, accelerates and deforms 
fluid blobs in a manner that prevents them from being compressed. The pressure 
is given by the solution of a Poisson equation, implying that it is a quantity dom- 
inated by the large scales of the velocity field, but the source term is quadratic 
in the velocity gradient and leads to non- Gaussian statistics of the pressure field. 
Recent studies have shown that the pdf of the pressure field is highly skewed and 
has a long tail for negative fluctuations. These fluctuations correspond to intense 
vortices with radii of the order of the Kolmogorov length (Holtzer & Siggia 1993, 
Pumir 1994). 

The internal dynamics of a turbulent fluid motion can be extracted in a frame 
moving with the fluid (the Lagrangian point of view). At moderate to high Reynolds 
number, the pressure gradient rather than the viscous stress is responsible for the 
deformation and acceleration of fluid regions larger than the dissipation scale. For 
example the time scale of small scales of turbulent motion can be inferred from the 
curvature of the Lagrangian velocity auto- correlation obtained from the expansion 

Ri{i,s) = (v(<) • v(s)) = Co - - s) 2 + •••, t>s, (1.1) 
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C 0 =2 f E(k, s)dk, 

Jo 

C 2 = — ((Vp(x,s)) 2 > = - jf°° k 2 P{k,s)dk = -j WM)0(-M)} 


( 1 . 2 ) 

k2 dk, (1.3) 

D(k) = kik t I dpdq u,(p, a)u/(q, s), (1.4) 

J J p+q=k 

for small t — s. Here for simplicity we have neglected the viscous term and have 
taken p — 1. The power spectrum of the pressure P(k,t) gives the pressure variance 
as 


(. p 2 (x,<))=/ P{k,t)dk. 

Jo 


Eqs. (1.3) and (1.4) imply that the spectrum of the pressure gradient depends on a 
fourth-order moment of the velocity field, and roughly speaking the largest contri- 
bution to the pressure-gradient variance comes from wavenumbers lower than the 
peak wavenumber of k 2 E(k ) when Reynolds number is large. This can be compared 
with the variance of the pressure itself, which is dominated by wavenumbers near 
the peak of the energy spectrum E(k), and with the variance of the dissipation e, 
which is dominated by wavenumbers near the peak of k 2 E(k). 

It is well known that the small scales of turbulent motion are intermittent (we 
shall use this term herein to mean simply a departure from Gaussian statistics) and 
fourth-order moments have contributions from the cumulant part. An example is 
seen in Gotoh et al (1993) where the initial curvature | is larger than |C^| 
computed from a Gaussian velocity field 


C? = ~<(Vp(x,5)) 2 ) = - / k 2 P G (k,s)dk 

Jo 

rOO r OO 

= — j dk j dq kq J E(k, s)E(q, s), 


(1.5) 


J(x) = {(a 2 - l) 2 log _ 2a + ^a 3 }/(2a 4 ), a = . (1.6) 

This suggests that turbulence has a faster decay of Rl{t) and a smaller turbulent 
diffusivity (as given by the time integral of Ri{r)) than predicted by the Gaussian 
theory. In other words, it suggests that the dynamics of the Navier-Stokes equation 
causes a fluid blob to forget its past history more effectively than would convection 
by a Gaussian velocity field. Therefore intermittence effects appear in the dynamics 
of a fluid particle through the pressure gradient and ore expected to be different 
from those of the viscous stress, but we do not know how fast intermittence in the 
pressure increases with wavenumber, how it differs from the intermittence of the 
dissipation, what its time scale depends on, and so on. 

There have been many studies of the fluctuations of pressure and its gradient in 
turbulence (Monin &; Yaglom 1975, and Nelkin 1994) but few fundamental studies 
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FIGURE 1. The pdf of pressure, (a) turbulent velocity field. : R\ = 172, 

: R\ = 96, : R\ = 63, — ■ — : R\ = 38, o: Gaussian, (b) Gaussian 

velocity field with the same energy spectrum for R\ = 172. : turbulent 

velocity field, H : Gaussian velocity field, o: Gaussian. 
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of pressure intermittence and its effects on the flow dynamics. This is partly because 
experimental measurement of the pressure spectrum is difficult (Uberoi 1953 and 
George et al 1984). Now that extensive numerical data are available, we can examine 
the statistics of the pressure field systematically, albeit at rather low Reynolds 
numbers. Here we present some measurements taken from the DNS data base 
available during the summer program and examine the statistics of the pressure 
field. 

2. Numerical simulation 

The turbulent flow fields used here (Jimenez et al 1993) are all homogeneous and 
isotropic and are held stationary by forcing at low wavenumber. They are arranged 
into four groups according to the Reynolds number R \ , with statistical quantities 
computed as averages over two fields at well separated times for R\ = 172 and 
over three fields for R\ = 96, 63, 38. The numerical grid sizes were N = 256 3 for 
R\ = 172,96, N — 128 3 for R\ = 63, and N = 64 3 for R\ = 38. The Lagrangian 
autocorrelations of velocity and pressure gradient were computed by the passive- 
vector method (Kaneda &; Gotoh 1991, Gotoh et al 1993). The pressure field p(k,f) 
was computed using de-aliased spectral methods. 

3. Results 


3.1. One-point one-time statistics 

The one-point pdf of the pressure in Fig. la is skewed as reported by Pumir 
(1994). For negative pressure fluctuations the asymptotic form of the pdf tends 
to be nearly exponential P(p) oc exp (-a\p/<j p \ a ) , where a is a non-dimensional 
constant and the exponent a is slightly less than one, which is consistent with low 
pressure in the core regions of intense vortex filaments. The tail of the pdf extends 
towards negative values with large amplitude as R\ increases but a is independent 
of R\. For positive fluctuations the pdf is close to Gaussian and is insensitive to 
R\ (Pumir 1994). Furthermore the pdf Pg(p) of the pressure field computed from 
a Gaussian velocity field having the same energy spectrum has the same behavior 
(Fig. lb) for the negative fluctuations as predicted by Holtzer & Siggia (1993). 

The pdf of one component of the pressure gradient is shown in Fig. 2a for different 
Reynolds numbers. The pdf is symmetric and its tails become wider as R\ increases. 
Fig. 2b shows that the pdf of the pressure gradient is isotropic and that it differs 
markedly from that computed from a Gaussian velocity field, unlike the case of 
the pressure itself. The asymptotic form of the tails of these pdf’s appears to be 
a stretched exponential P(p } i ) ex exp b \p,i/&p ti with 0 < ft < 1, where b is 
a non-dimensional constant, while that of its Gaussian counterpart is exponential 
as expected. This means that the pressure-gradient field is very intermittent and 
is quite different from the pressure gradient in a Gaussian velocity field. Holtzer Sz 
Siggia (1993) studied the pdf of the pressure gradient and suggested ft ~ 1/2. The 
pdfs plotted against y/p t yf<Jp ~ (figure not shown) indicate that /? is close to 1/2 for 
large amplitudes and becomes smaller as R\ increases. 
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FIGURE 2. The pdf of pressure gradient, (a) variation with Reynolds number 

for a single component. : R\ = 172, : R\ = 96, : R\ = 63, 

: R\ = 38, o: Gaussian, (b) Comparison of turbulent and Gaussian values 

for all components at R\ = 172. : dp/dx, : dp/dy, : dp/dz for 

turbulent velocity field; H : dp/dx, □ : dp/dy, x : dp/dz Gaussian 

velocity field; o: Gaussian. 




194 


T. Gotok & R. S. Rogallo 




FIGURE 3. Variation with Reynolds number of velocity moments, (a) Various 
moments of turbulent velocity field. 1: F p , 2: JFVp, 3: F c , 4: 5: ttl , 6: 

F ttl . (b) Comparison with values from the corresponding Gaussian velocity fields. 
1: F p (DNS), 2: F p (Gaussian), 3: F p (George et al 1984), 4: Fv P (DNS), 5: Fx? P 
(Gaussian), 6; FVp (George et al 1984). 


The variation with R\ of normalized moments of the pressure and pressure gra- 
dient are shown in Fig. 3a. Included there for comparison are other normalized 
fourth-order moments of velocity: 


(p ! > 


F " aw 


Fv P = 


((Vp) 2 ) 

p/ij,- 1 / 2 ’ 


(f!> 

<<> 

where e = (c) 


F( — , 1 2 > 


F„ = 


4 ) 


— 


((dui/dsi) 4 ) 




((m) 4 ) 


( 3 - 1 ) 


(3.2) 


< W 2) 2 ’ ^ lUl {(^u^axj) 2 ) 2 ’ ttl <( Ul ) 2 ) 2 ’ 

The normalized variance of the pressure F p is insensitive to R\ 
(Batchelor 1951 and Pumir 1994). On the other hand the normalized variance of 
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k/kd 

FIGURE 4. Velocity-gradient and pressure- gradient power spectra in Kolmogorov 

units. : R\ = 172, : R\ = 96, : R\ = 63, : 

R\ = 38. Curves without symbols: e~ 3 ^u 1 ^k 2 E(k). Curves with symbols: 
/4jfc 2 P(Jfc). 

the pressure gradient Fv ? increases rapidly with R\ and is roughly proportional to 
R 1 ^ 2 for the range of Reynolds numbers studied here. This tendency is in agree- 
ment with Yeung Sz Pope (1989), where R\ was below 93. On the other hand F € 
and Fu increase slowly with R\ in agreement with the data of Kerr (1985). This 
difference is surprising because functions of the velocity gradient such as vorticity or 
du\/dx\ axe quantities representing small scales of turbulent motion (characterized 
by wavenumbers near the maximum of k 2 E(k)) and are believed to be more inter- 
mittent than the pressure-gradient field whose characteristic wavenumber is lower 
than that of k 2 E{k) when the Reynolds number is high (see Fig. 4). 

When the universal equilibrium form of the energy spectrum (Kolmogorov 1941) 
is used with the Gaussian-velocity approximation (1.5) we obtain 

E(k) = i l W 4 f(k/k d ), 

((Vp(x)) 2 > = (33) 

Fv P = 2 f dx j dy xy J (^) f(x)f(y), 

Kaneda (1993) analyzed (3.3) and found that the largest contribution to the integral 
comes from the region y < k c /kd < x, where k c is the wavenumber above which 
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E(x ) oc exp(— cr) and c is a non-dimensional constant, which means that the integral 
Fv P is governed by the energy spectrum up to the dissipation range. Batchelor 
(1951) estimated Fv P ~ 3.9 independent of R\. George et a l (1984) used (3.3) with 
an empirical energy spectrum for E(k ) and found 


Fv P = 



(3.4) 


with flp = 3.7 and 7 P = 62.7. The constant f3 p is universal while 7 P depends on the 
macro-scale of the turbulence. 

Fig. 3b compares the variances of the pressure and its gradient with those com- 
puted from Gaussian velocity fields. The values computed from the Gaussian fields 
are close to those found from (3.4) using the values of George et al (1984) but are 
significantly lower than the DNS values. The R\ dependence of the turbulent pres- 
sure gradient is stronger than that of its Gaussian counterpart, and it seems likely 
that any theoretical explanation must take into account the non- Gaussian statistics. 

3.2. Two-point one-time statistics 

Eqs. (1.3) and (1.4) show that the power spectrum of the pressure gradient is a 
function of a fourth-order moment of the velocity field. The Gaussian approximation 
for the velocity field leads to (1.5), which implies that the peak of k 2 P(k ) occurs at 
a lower wavenumber than that of k 2 E(k). Comparison of these spectra, plotted in 
Kolmogorov units in Fig. 4, confirms that at low R\ both spectra peak at nearly 
the same wavenumber, but as R\ increases the peak of k 2 P(k) moves to lower 
wavenumber. Collapse of the enstrophy spectra is excellent but that of the pressure 
gradient is not. This implies that the Kolmogorov scaling is not appropriate for 
the pressure field (Moin & Yaglom 1975). The cumulative contributions of these 
spectra to the total enstrophy and pressure-gradient variance 


Qi(k) 


Jo* ^Ejg)dq s /o* ^Pjo)dq 

/*”•“ q*E(q)dq * /*— q*P{q)dq 


(3.5a, b ) 


for the four R\ s reach 80% at k p /kd « 0.35 for the pressure gradient and at 
ku/kd & 0,5 for total enstrophy, but their scale separation k p /ku is small (Batchelor 
1951). 

One way to measure the variation of pressure intermittence across the spectrum 
is to compare the pressure power spectrum P(k) of the turbulence with that Pg(&) 
computed from a Gaussian velocity field (Hudong et al 1987). Fig. 5a presents such 
a comparison for R\ = 172. The pressure spectrum P(k) is larger than Pg(^) for 
all wavenumbers beyond the forced range, meaning that the cumulant contribution 
Pc(k ), 

P(k) = P G (k) + Pc(k\ (3.6) 

is positive. This suggests that the inertial-range constant K p for the pressure spec- 
trum P(k) = A'pe 4 / 3 ^" 7 ^ 3 is larger than that predicted from the Gaussian velocity 
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Figure 5. Comparison of pressure power spectra in turbulence with those com- 
puted from Gaussian velocity fields, (a) at R\ = 172. : P(k), : Pc(k). 

(b) variation of their (scaled) ratio R^ 1 ^ 2 K\(k) with Reynolds number. : 

R x = 172, : R x = 96, : R x = 63, : R x = 38. 
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field (Nelkin &; Tabor 1990, Fung et al 1992, and the article by Pullin & Rogallo in 
this proceedings). We consider the spectrum ratio 


K\(k) = 


m 

Paik) 


(3.7) 


as a measure of the variation of pressure (and pressure gradient) intermittence across 
the spectrum. A R x ' dependence of ivi(fc) is suggested by the good collapse of 

_ 1/2 w 9 B 9 4 ^ 

i? A ' A"i(A;) shown in Fig. 5b. This implies that 

Ki(k) = 1 + cx -iZi /2 log (A) , */*„ < 1/2. (3.8) 

The growth with wavenumber of the intermittence of second-order moments of the 
pressure field is very slow (it appears to be logarithmic, but the scale range is not 

sufficient to preclude a weak algebraic dependence) but the growth with Reynolds 

1 /2 • • » 
number as RT * faster than that of other small-scale quantities. At R\ = 172 


((Vp) 2 ) 

<(Vp)% 


2.58 


(3.9) 


and the initial curvature of the Lagrangian velocity autocorrelation is larger than 
that of the Gaussian field (Gotoh et al 1993). Navier-Stokes dynamics decorrelates 
the velocity of fluid particles faster than convection by a Gaussian velocity field 
would. 

The variation of the dissipation intermittence across the spectrum is measured in 
a similar way. Fig. 6 compares the dissipation power spectra of turbulent-velocity 
fields with those of the corresponding Gaussian- velocity fields. Here we define the 
power spectrum of the dissipation and its intermittence (ratio of turbulent to Gaus- 
sian spectra) as 


(f 2 (x,0) = f 
Jo 


E € (k,t ) dk , 


(3.10) 


K 2 (k) = 


EJJc) 

E ( ,a(kY 


(3.11) 


The intermittence within the dissipation range is independent of Reynolds number. 
It is interesting to note that the maxima of J\Ti ( Ar) and minima of K 2 (k) both occur 
at k/kd ^ 1/2 and that the diameter of intense vortex tubes is also of order 1 jk^ 
(Jimenez et al 1993). The strong variation of the dissipation intermittence across 
the spectrum is opposite to that of the pressure field due to the modulation of the 
energy spectrum by the small-scale structures (Hudong et al 1987). For example 
when strong singularities of the velocity field with support size 1 fkd are placed 
periodically at a fixed separation l 1/Arj, the power spectrum of the dissipation 
field has an excitation at k ~ l/l <C kj as well as at k ~ kj, while a Gaussian 
velocity field with a compact spectrum centered on kj simply generates a compact 
dissipation spectrum centered on k j. 
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FIGURE 7. Lagrangian correlations of velocity and pressure gradient, (a) R\ = 63. 
(b) R x = 96. : R$ p (t, s), : R$(t, s). 



Figure 8. Kolmogorov time-scaling of the Lagrangian pressure-gradient correla- 
tion. : R\ = 96, : Rx = 63. 
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FIGURE 9. The decay of Lagrangian and Eulerian pressure-gradient correlations 

within octave bands at R\ = 63. : Lagrangian, : Eulerian. The 

uppermost curves correspond to the lowest octaves and the correlations decrease 
monotonically for higher octaves. 


3.3 . Two-time statistics 

The normalized Lagrangian auto-correlations of velocity and pressure gradient 
are defined in a stationary flow as 


_ (v(x,*10 • v(x,*|a)) 
Kt ) (Kx,3|3)| 2 ) 

(F(x,*|0-F(x t <ls)) 

’ } (|F(M*)l a > 


•Rv P (* 


for t > s, 
for t > s y 


(3.12a) 

(3.126) 


where v(x, t\s) is the velocity at time s of the fluid particle whose space-time tra- 
jectory passes through (x,£), and F(x,t|s) is the pressure gradient acting on that 
particle at time s. Note that v(x,t|t) = u(x,t) and F(x,t|t) = Vp(x,t). 

Fig. 7 compares the Lagrangian auto-correlations J?^(t,s) and i?y p (f,s) at R\ = 
63 and 96. The correlations of the pressure gradient decay much faster than those 
of the velocity and have negative correlation at later times. In Fig. 8 R^ p (t y s) is 
replotted with time scaled in Kolmogorov units. Collapse of the curves is excellent 
up to the time at which they change sign. 

We define octave-band Lagrangian correlation spectra for the pressure gradient 
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FIGURE 10. Time-scaled narrow-band Eulerian correlations of velocity and pres- 
sure gradient for wavenumber-bands k = 2 + 2n,n = 0, ■■■29 at R\ = 63. 
(a) R*(k,t, S ). (b) R% p (k,t,s). 


as 


3 ) — 


£(F(M<)-F(-k,<| 5 )) 

oct 

52 |F(k,t|t )| 2 52 |F(k,<|s )| 2 

J oct oct 


1/2 ’ 


(3.13) 


where 53 = 53 *= \ + and F(k, t\s) is the Fourier transform of F(x, t\s ) with respect 

oct 

to the Lagrangian coordinates at labeling time t (Kraichnan 1966 and Gotoh et al 
1993). The octave-band Eulerian correlations are defined in a similar way. 

The comparison of these spectra in Fig. 9 indicates that at high wavenumbers the 
Lagrangian correlation decays slower than the Eulerian, but that at low wavenum- 
bers the decay rates are roughly equal. The Eulerian decorrelation of the small 
scales is caused primarily by their sweeping by larger scales while their Lagrangian 
decorrelation is due to deformation by larger scales. The sweeping effect on the 
Eulerian correlation is clearly demonstrated in Fig. 10 where narrow-band Eulerian 
correlation spectra, defined similarly to (3.13) with 52 0 ct replaced by 53aAef/> BTe 
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FIGURE 11. Time-scaled narrow-band Lagrangian correlations of velocity and 
pressure gradient for wavenumber- bands k = 2 + 2 n y n = 0, ■••29 at R\ = 63. 
(a) R„(k,t,s). (b) R^ p (k y t y s). The low-wavenumber bands display pronounced 
negative values. 


plotted against the normalized time kUo(t — s). On the other hand, in Fig. 11, 
the narrow-band Lagrangian correlation spectra collapse when plotted against the 
normalized time r)(k)(t — s), where 


V (k) = 



(3.14) 


is the rms strain rate of eddies larger than 1/k (the inverse of the turnover time 
for an eddy of that size), which is obtained from (1.5) by using the expansion 
of J(x) ~ Ijifx for small x (Gotoh et al 1993). The success of these scalings is 
confirmed by a comparison of Figs. 10b and lib with the unsealed data in Fig. 9. 


4. Conclusions 

Pressure statistics have been examined within a DNS of stationary isotropic tur- 
bulence. The one-point pdf of the pressure has an exponential tail for negative 
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fluctuations while that of the pressure gradient is symmetric and of stretched ex- 
ponential form. The normalized variance of the pressure is insensitive to R\ but 
the normalized variance of the pressure gradient increases as iZ A ' * The variances 
of both pressure and its gradient in the DNS are higher than those computed from 
a Gaussian field having the same energy spectrum. 

The growth of intermittence across the pressure spectrum, characterized by the 
ratio of power spectra P(k)/ Po(k), was found to be proportional to — R x f log (k/kd) 
for k/kd <1/2. The Lagrangian pressure-gradient correlation has a time scale much 
shorter than that of the velocity and its narrow-band spectra change sign at the 
normalized time 77 ( Ar ) ( ^ — s) ~ 1. This is consistent with the notion of a fluid particle 
being accelerated toward the center of a vortex of size 1/k by the pressure gradient 
acting on that particle. 

Previous theoretical analyses of these observations are rather scarce (Pullin 1994 
and Nelkin 1994) due to the difficulty of dealing with the nonlocality in physical 
space of the Poisson equation for the pressure field. Most studies so far are based 
on the Gaussian velocity approximation, but intermittence effects in the pressure 
field are clearly evident in its one-point statistics a s well as its spectrum. The 
construction of a theory that is able to explain these observations is a real challenge. 
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